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ABSTRACT 

We have used the Berkeley-Illinois-Maryland Array (BIMA) to make deep N2H+ 1-0 maps 
of the pre-stellar core LI 83, in order to study the spatial and kinematic substructure within 
the densest part of the core. Three spatially and kinematically distinct clumps are detected, 
which we label L183-N1, L183-N2 and L183-N3. L183-N2 is approximately coincident with 
the submillimetre dust peak and lies at the systemic velocity of L183. Thus we conclude that 
L183-N2 is the central dense core of L183. L183-N1 and 3 are newly-discovered fragments 
of L183, which are marked by velocity gradients that are parallel to, but far stronger than, the 
velocity gradient of LI 83 as a whole, as detected in previous single-dish data. Furthermore, 
the ratio of the large-scale and small-scale velocity gradients, and the ratio of their respective 
size-scales, are consistent with the conservation of angular momentum for a rotating, collapsing 
core undergoing spin-up. The inferred axis of rotation is parallel to the magnetic field direction. 



'Department of Astronomy, University of Illinois at Urbana-Champaign, 1002 West Green Street, Urbana, II 61801 



-2- 



which is offset from its long axis, as we have seen in other pre-stellar cores. Therefore, we 
propose that we have detected a fragmenting, collapsing, filamentary, pre-stellar core, rotating 
about its B-field, which is spinning up as it collapses. It will presumably go on to form a 
multiple protostellar system. 

Subject headings: ISM: Individual: Alphanumeric: L183 — Stars: Formation — Radio Lines: 
ISM 



1. Introduction 



The star formation process begins within molecular clouds, where the internal pressures (magnetic, 
turbulent, thermal, or otherwise), suppo rting dense regions against co llapse, lose out to gravity. The densest 
regions are known as pre-stellar cores dWard-Thompson et all 1 1994b . which are the gravitationally bound 
cores of clouds that aie beli eved to be undergoing star fo rmation, but as yet show no signature of an em- 
bedded protostellar source (| Ward-Thompson et all 120070 . Recent studies have revealed the presence of 
extrem ely faint protostellar objects inside a srn all number of cores that were previously believed to be pre- 
stellar (|Crapsi et al.ll2005bl : iBourke et al.lbOOa') . These appear to be very young protostars embedded deep 
within these cores. 

However, the manner of the collaps e from a pre-stell a r core to a core that contains a protostar is still a 
matter of some debate - for a review, see IWard-ThompsonI (120021) . Gravitational collapse, rotation, turbulent 
motions and fragmentation may all play a role. Detailed study of the densest pre-stellar cores that do not 
contain a protostar is necessary to understand the evolution of pre-stellar cores (IJessop & Ward-Thompson 

booik 

The centrally condensed, high column density pre-stellar core LI 83, which is also known as L134N, is 
part of a loose group of high latitud e (b ^37°) de nse molecular cores associated with the surface of the local 
bubble at a distance of 1 10 =fc 10 pc (IFrancdl 19891). It has be e n extensively studied with a number of different 



spectroscopic spec ies (|Lee et all 1200 ll : iDickens et al.ll2000l : lPagani et al.ll2005l: iReqi 
Pagani et al .I2OO7I) a nd in the far-infrared and submillimetre cont i nuum Jjuvela et al.l 



Reauena-Torres et al 



2OO2I : Lehtinen et al]l2003l : IPagani et al.ll2003l : 



Kirk et al 



2005; 



Kauff mann et al. 



2007 



2OO2I : I Ward-Thompson et al 



20081) . It also appears to 



be more chemically evolved than other pre-stellar cores ( Dickens et al.„200(]i) and displays a number of sig- 



nature s common to cores that are believed to be close to forming a protostar (|Lee et al.ll200ll : ICrapsi et al. 



2005 al) . Thus, LI83 is a good candidate in which to study pre-stellar evolution. 



In this paper we present new N2H+ observations performed with the BIMA interferometer. N2H+ is 
one of the preferred tracers for cold dense molecular material as it is well correlated with submillimetre 
dust emission and is less prone to gas phase depletion than other molec ular species within pre-stellar cores 
(|Jessop & Ward-Thompsonll2001l : lBergin et al iTafalla et"aDl2004l), and it has been widely used to sur- 
vey pre-stellar cores in ge neral and L183 in particular ( Dickens et al. 2000; Caselli et al. 2002; Crapsi et aP 



2005al : IPagani et al.ll2007l ). 
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In Section |2] we describe the configuration of the BIMA array and in Section [3] we show channel and 
integrated maps of N2H"'" intensity. In Section |4] we show the results of fitting various line parameters, and 
we discuss their implications in Section [51 



2. Observations 



BIMaB was used t o map the 3.22 mm (93 GHz) J = 1 ^ transition of N2H+ Jcaselli et al. 



1995 : Daniel et al. 20o d) towards the peak of the L183 submillimetre dust emissioiH. BIMA was a 10- 



element millimeter -wave interferometer operated at the Hat Creek Radio Observatory in Northern California 



dWelch etal.ll 19961) . 



BIMA was used in its compact C-array (~2-27 kA) and ultra-compact D-array (~2-9 kA) configura- 
tions. Table [T] shows the individual observation dates, durations and configurations. A total duration of 50 
hours (~6 tracks) was spent in the C-Array configuration and a total duration of 18 hours (~2 tracks) was 
spent in the D-Array configuration. 

Our observations are only sensitive to angular scales between the FWHM of the primary beam (given 
by the effective diameter of an individual dish) and that of the synthesized beam (given by the uv sampling 
and weighting). The individual antennae have a physical diameter of 6.2 m and an effective primary beam 
FWHM of 1.7 arcmin at 3.22 mm. The synthesized-beam FWHM for the combined observations was 13 x 10 
arcsec with a semi-major axis position angle of —14°. The beam was elongated due to the declination of the 
source. Therefore, these observations are sensitive to structure between 10 and 100 arcsec. This corresponds 
to a spatial scale of 0.005-0.05 pc at 110 pc. 

The BIMA correlator was set to mode 5 to give 512 channels across a bandwidth of 6.25 MHz. The 
channel width was 12.2 kHz (0.039 km s^^). The N2H+ lineset was placed in the first sideband while 
the 96.4 MHz C^^S 2-1 transition was in the opposite sideband, but no signal was detected in the second 
sideband. Phase and amplitude calibration were conducted against the radio source 1549-1-026. The assumed 
flux of 1549-1-026 at our wavelength was 1.5 Jy (BIMA Calibrator Lis@). The data were reduced using the 
MIRIAD software package and the visibilities were Fourier tran sformed into the image plane using natural 



weighting in order to maximize the sensitivity (ISault et al.lll995h . The theoretical per channel rms was 0.121 
K. 



'Operated by the Berkeley-IUinois-Maryland Association with NSF and University Funding. The BIMA facilities have since 
been merged with the OVRO array to form the present CARMA array. 

^a2ooo=15'^54™08.8^ (52ooo -~02°52'38" [Kirk et allboosh . 

'http://bima.astro.umd.edu/cgi-bin/swflux_plot.pl?source=1549%2B026 
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3. Results 



The total intensity of emission from the seven hyperfine components of N2H+ towards the LI 83 core 
is shown in Figure [T] integrated over velocity. The first solid contour is at 2.46 Jy/Beam, which is also the 
separation between subsequent contours. For this dataset the scaling between Jy/Beam and Kelvin is 1.07 
K/(Jy/Beam). 

The integrated emission shown in Figure [T] displays the characteristic north-south L183 filament as 
seen in submillirnetre dust emission, infra-red extinct ion, and single dish line maps of N2H"'" 1-0 (e.g. 
Pagani et al.ll2003l : iKirk et al.ll2005l : ICrapsi et al.ll2005ah . The integrated N2H+ intensity shows three bright 
peaks (north, south and west). For reference we will refer to these features as LI 83 N2H+ (L183-N) 1, 2, 
and 3 in order of their peak intensity. Their positions are listed in columns 2 and 3 of Table |2] and are shown 
in Figure [Hand in subsequent figures as labeled open stars. 

To examine the kinematic structure of LI 83 we extracted channel maps of N2H+ intensity across the 
101-012 component as 0.039 kms^^ slices between 2.15 and 2.72 kms~^. This component is shown as 
it is isolated in velocity from the other components and thus is free from blending between components. 
These channel maps are shown in Figure [2] As the velocity increases, the emission first appears in the 
northeast, before moving to the south, and finishing in the northwest. These components roughly correspond 
to the three emission peaks seen in Figure [T] Thus, L183-N1 and L183-N3 are separated spatially and 
kinematically from each other, but are both connected kinematically to L183-N2. 



The central section of the SCUBA 850-/um map from ' Kirk et al.l(|2005h is shown as a greyscale image in 
Figure[3] The SCUBA continuum emission peaks close to L183-N2 and has a northward extension betwee n 
L183-N1 and 3. This extension is seen to continue on a larger scale in the SCUBA image dKirk et allbooi) . 
The L183-N2 peak is not exactly coincident with the SCUBA 850-/im emission peak. The peak in the 
continuum emission is in fact coincident with the N2H"'" emission at 2.44-2.48 km s^^. 



4. Spectral Line Analysis 



The data shown in Figure HJuse a pixel scale of 1.87 arcsec per pixel in order to show a smooth intensity 
distribution. For the analysis of the spectral lines the data were spatially averaged onto a 5.6 arcsec pixel 
grid in order to produce pixels that are approximately Nyquist sampled and to generate a factor of 3 increase 
in signal-to-noise ratio. The spectrum for each pixel that had a 2a or greater feature at ~2.5 kms^^ (the 
approximate location of the 101-012 hyperfine component) was th en extracted and piped to the hyperfin e 



structure (hfs) routine from the CLASS data reduction package^ ( Hily-Blanl 2000 : Buisson et al. 



2005). 



The hfs routine was used to simultaneously fit all seven hyperfine components giving best fits for the radial 
velocity Visr, the linewidth /S.V, the total optical depth r and the product of the main beam temperature 
and optical depth {Tmb x 't)- For the fitting it was assumed that all hyperfine components had the same 



"•class is part of the GILDAS software suite available from |http: / /www, iram. f r/IRAMFR/GILDAS| 
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Table 1 : Observation dates, duration and array configuration for the data used in this paper. 



Date 


Duration 


Configuration 


2003 March 23 


8 hrs 


C-Array 


2003 March 24 


8 hrs 


C-Array 


2003 May 03 


5 hrs 


C-Array 


2003 May 07 


8 hrs 


C-Array 


2003 May 08 


4 hrs 


C-Array 


2003 May 15 


5.5 hrs 


C-Array 


2003 July 16 


7 hrs 


D-Array 


2003 July 17 


7 hrs 


D-Array 


2004 April 30 


3 hrs 


C-Array 


2004 May 09 


8 hrs 


C-Array 


2004 May 14 


4 hrs 


D-Array 



Table 2: Table of positions and hyperfine fit parameters towards the three peaks in the integrated N2H+ 
map (Figure [D. Column 1 lists the source name. Columns 2 and 3 list the central position of each source. 
Columns 4 to 7 list the CLASS hfs parameters - radial velocity (Visr), linewidth (AV), total optical depth 
summed across all seven hyperfine components (r), and the main beam temperature (Tmb) - for a fit to the 
spectra at the positions listed in columns 2 and 3. Column 8 lists the non-thermal linewidth (AV/vt) based 
on the assumption that = 7 K (see Section |42l ). Column 9 lists the rms of the radial velocity (Vrms) 
computed in a 5 by 5 pixel box. Column 10 lists the excitation temperature (Tex)- Columns 11 and 12 list 
the magnitude and direction, and their standard deviations, of the local velocity gradient (see Section [5]l. 



Source 


R.A. 


Dec. 


V,,r 


AV 


T 


Tmb 


AVnt 




T„ 


S 


PAg 




(2000) 


(2000) 


[kms-1] 






[K] 


[kms-i] 




IK] 


[km s^^ pc^^] 


[°] 


L183-N1 


15''54™09.5^ 


-02°52'15" 


2.256±0.003 


0.123±0.006 


0.37±0.08 


1.2±0.4 


0.064 


0.054 


6.5±2.0 


6.5±5.8 


-74±23 


L183-N2 


15h 540199 js 


-02°52'50" 


2.415±0.004 


0.127±0.008 


2.1±0.6 


0.5±0.2 


0.071 


0.025 


3.3±1.8 


1.0±2.4 




L183-N3 


15^54^07.3^ 


-02°52'31" 


2.588±0.005 


0.108±0.010 


0.36±0.18 


0.63±0.45 


0.025 


0.090 


4.8±3.8 


8.9±8.0 


-57±60 
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15^54"'10^ 08^ 06^ 

Right Ascension (2000) 



Fig. 1. — A map of N2H+ (1-0) line emission integrated across the seven hyperfine lines towards the 
L183 core. The greyscale shows the full range of emission, the contours start at, and have intervals of 2.46 
Jy/Beam. The dashed contours are negative, they start at, and have intervals of —2.46 Jy/Beam. The hatched 
ellipse shows the size of the synthesized beam. The three sources, L183-N1, 2 & 3, are marked with stars 
and are numbered. The contour spacing is 1.5 times the off-source rms of 1.64 Jy/Beam. 
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15''54°10' OB^ 06^ 15''54"10' 08° 06' 15''54"l0' 08" 06' 15''54°10' OB' 06' 

Right Ascension (SOOO) Right Ascension (SOOO) Right Ascension (2000) Right Ascension (8000) 



Fig. 2. — Channel map across 2.15-2.72 km s^^ of the N2H+ 101-012 hyperfine component towards L183. 
The channel width is 0.039 kms~^. The velocity of each channel is shown in the top-left corner of each 
panel. The contours are in increments of 0.156 Jy/Beam (0.169 K) with the dashed contours mirroring neg- 
ative values of the positive solid contours. The stars are as in FigurelH but only the stars in the 2.13 km s~^ 
panel are labeled. The contour spacing is 1.5 times the off-source rms of 0.104 Jy/Beam. 
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Fig. 3. — Inte grated NgH "* " cont ours from Figure [T] overlaid on a greyscale image of SCUBA 850-^m dust 
emission from iKirk et all ( 2005 ). The dotted-contour is at 0.16 Jy/Beam and traces the central part of the 
SCUBA emission. The stars and labels are as in Figured] 
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linewidth and excitation temperature. 

Table |2] lists the fit parameters at the positions of the three peaks identified from the map of integrated 
emission. The spectra for these positions are shown in Figure |4l We note that the errors on the fitting of 
Visr and Ay are narrower than the channel width (0.039 km s~^), due to multiple channels over the line that 
contribute to the fitting. 



4.1. Kinematics 



Figure [5] shows the variation of the local standard of rest velocity (Visr) across the L183 molecular 
core. It shows the same triple-lobed structure as shown in Figure [2l but now we see that the most extreme 
velocities are towards the outside of each of the three lobes. The velocities at the position of each of the 
three L183-N sources are Usted in column 4 of Table|2l There is a 0.33 km s^^ separation between L183-N1 
andL183-N3. 

There have been previous observations of N2H" '" towards L183 usin g the single-dish IRAM 30-meter 
telescope. Its beamwidth at 3.22 mm is ~30 arcsec. IPagani et al.l (|2007h measured a value of V^r = 2.367 
kms~^ towards RA(2000) = 15^5 4"^08.5^ Dec.(200q ) = -02°52'48". This is in good agreement with 
the values we measure in that area. Icrapsi et al. ( 2005a ) measured a value of Visr = 2.413 kms~^ towards 
RA(2000) = 15'^54"^08.4^ Dec.(2000) = -02°52'23". This again is similar to the values we measure 
exactly at that point. The literature value for the Visr of the whole core of L183 (2.415 km s"^) matches 
closely that of L183-N2. Hence we take 2.415 km s~^ as the systemic velocity of the LI83 core and we see 
that the velocities of L183-N1 and L183-N3 are velocity structures overlaid on this, of up to ±0.17 kms~^. 

To examine the velocity gradient across LI 83 we averaged the gradient vectors between a single pixel 
and its nearest neighbors with the condition that each pixel must have at least seven neighbors with valid 
measurements of Vsr- A vector map of the local velocity gradient is overlaid on the map of Vsr in Figured 

Figure 5] also shows a blac k-line indicating the direction of the plane-of-the-sky magnetic field as 



observed by lCrutcher et al.l(l2004b . This is approximately orthogonal to the large-scale velocity gradient (see 
Sect ion HI). The magneti c field may therefore play an important role in the evolution of this pre-stellar core 
(c.f. [Holland et al.lll996|) . We discuss further the relevance of the velocity gradient to previous observations 
in Section [5] 



4.2. Linewidtlis 



Figure [6] shows the remainder of the results from the hfs fitting. Panel (a) shows the linewidth vari- 
ation across LI 83. The linewidths measured from the hfs fitting are narrow. The mean linewidth in Ta- 
ble [2] is only 0.119 kms~^ and this is almost within the la dispersion of each measu rement. The htera- 
ture values for N^H"^ linewidths in LI 83 on peak are in the range 0.195-0.25 kms"^ dCaselli et al.ll2002l : 



Crapsi et al.l l2005al : IPagani et al.ll2007h . which in turn are on the narrow side of the average linewidth of 
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L183-N3 




Velocity (km s ') 



Fig. 4. — The N2H+ spectrum (gray line) observed observed towards the L183-N1, 2, and 3 positions 
averaged over a 5.6 arcsec box. The black line is the CLASS hyperfine fit to each position. The parameters 
of each fit are listed in Table |2] 
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Fig. 5. — A map of radial velocity. The overlaid contours are the integrated contours from Figure [T] The 
range of the velocities is shown by the bar at the top of the plot. The three labeled stars show the positions 
listed in Table|2] The plot is annotated with a 0.01 pc lin ear scale bar and a blac k-line showing the orientation 
of the plane-of-the-sky magnetic field as measured by ICrutcher et all (12004) . The vectors show the local 
gradient of the Visr- The length of the vectors shows the size and direction of the velocity gradient with 1 
arcsec equal to 1 kms~^ pc~^. 
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0.250 ± 0.07 kms-i for the 31 starless cores surveyed by lCrapsi et all (l2005ah . 



A ridge of broader linewidths of about 0.14 kms~^occurs between L183-N2 and the northern part of 
the core containing L183-N1 and 3. The broadest linewidth along this ridge occurs between L183-N1 and 
L183-N3. This structure can be attributed to blending between the separate velocities of L183-N1 and 3. 



A typical observed spectral line can be decomposed into non-thermal AVnt and thermal AVr com- 
ponents where the total linewidth is AV"^ = AV^j, + AK^. The thermal line width is given by AVr = 
Cs a/8 In 2. The sound speed is given by = kTx / yi.rnH where Tk is the kinetic temperature, mn is the 
mass of a hydrogen atom, k is the Boltzmann constant, and /i is the mean particle mass in units of atomic 
mass number (for N2H"'" this is 29). 

From the above we can estimate the non-thermal component from knowledge of the local sound speed, 
but we first need an estimate of the kinetic temperature. We have previously estimated the dust temperature 
in LI 83, using a greybody fit to the extended infra-red and submilli metre emission (in a 15 0-arcsec diameter 
aperture) centered on the dust peak, and found a value of lOK dward-Thompson et al. Il2002h . If this were also 
the kinetic temperature of the gas, then the thermal linewidth would be 0.126 kms~^, and the non-thermal 
linewidth would be zero across the core. 

However, radiative transfer studies have shown that the extinction within pre-stellar cores can support 
very low temperatures at the center (a s low as ^ 7 K) with the tempe rature increasing radially outwards 
(lEvans et al.ll200ll : IZucconi et al.ll2001uStamatellos & WhitworthI 120031) . In this scenario the dust observa- 
tions give a mean temperatur e, but the temperature in the den sest central parts, as probed by the current 
observations, is usually lower (IStamatellos & Whitworthll2003 h. 



Another way t o calc ulate the temperature is to use the relative strength of different molecular transi- 



tions. IDaniel et al.l (l2007h and lPagani et all (120070 recently d id this for higher quantum number transitions 
of N2H"'" a nd its deuterated iso topologue. The IPagani et al.l (120071) model found an internal temperature 
of 7K. The iDaniel et al.l (l2007h model found an internal temperature of 8K within a central radius of 20 
arcsec. So we adopt the kinetic temperature of Tk=1 K. From this we calculate a thermal linewidth of 
AVt'=0.105 kms~^and the non-thermal linewidths listed in column 8 of Table|2] A kinetic temperature of 
Tk=^ K would give a AT/t=0.112 kms^^. 

Alternatively, we can calculate Vrms, the plane of the sky variation of Visr, under the assumption that it 
is produced by AVnt- This will be a measure of AVnt, but will not necessarily be equal to it. We calculate 
Vrms by calculating the root-mean-square of the radial velocity in a 5 x 5 pixel box centered on each of the 
three peaks (~ 30 x 30 arcsec). These values are listed in column 9 of Tabled 

The value of Vrms for L183-N2 is lower than for the other two sources. This might be expected from a 
study of Figure m in which it can be seen that the values of Visr appear fairly uniform across L183-N2. In 
the case of L183-N1 the mean values of AVnt and Vrms only differ by a few ms^^. The variation is slightly 
greater for L183-N2 and 3, but for all three sources the mean of AV^t and Vrms is about 0.05-0.06 kms^^. 

Particularly narrow linewidths are measured along the southern edge of the L183-N2 core. To the 
south of L183-N2 the total linewidth drops to 81 ± 6 ms^^. If this were purely thermal it would imply a 



temperature of ?V=4 K. However, we note that this is the region of highest optical depth (see next section), 
and so we may not be sampling the whole cloud along this line of sight. 



4.3. Optical Depths 



A map of the optical depth is shown in Figure [6J5. Column 6 of Table |2] lists the total optical depth 
summed over the seven hyperfine lines. The optical depth for any single line is equal to the total optical 
depth multiplied by th e fractional s trength of that line. As no single line contributes more than 25% to the 
total opacity, all lines (IDaniel et al.ll2006 ) will be optically thin for a value of r < 4. 



The two northern sources have total optical depths of order one-third while the southern source has an 
optical depth close r to two. The relatively low values of optical depth compared to earlier studies (20.3, 
Crapsi et all l2005al ') is due to the interferometric process resolving away the greater part of the extended 
material across this core. The only part that may be optically thick is the extreme southern edge. 



4.4. Temperatures 

The MIRIAD dataset was calibrated in Jy/Beam. For the CLASS hfs fitting a conversion factor of 1.07 
Jy/Beam/K was used to covert the extracted spectra into units of main beam temperature Tmb- Figure |6]: 
shows the distribution of fitted Tmb across the core. It is notable that the pattern of Tmb does correlate 
exactly with that of the integrated intensity. The intensity valley that separates L183-N3 from L183-N1 and 
2 is more pronounced. This corresponds to the northwestern edge of the emission in the 2.44 km s~^channel 
map and the southeastern edge of emission in the 2.60 kms^^channel map. 

The CLASS hfs routine returns the parameter x = Tmb x Thus the quoted Tmb error contains 
error contributions from the fit of x and r. These uncertainties also reflect the difficulty in fitting the height 
of the line when there are only a few channels across the peak. We can calculate the excitation temperature 
from the equation 

TMB = ^ = {Tex-Tbg){l-e-^) (1) 

where T^^ is the background temperature and is assumed to be equal to the cosmic microwave background 
(2.73 K). 

Figure|6]i shows a map of T^x- The values for Tmb and Tg^. at the L183-N peaks are listed in columns 7 
and 10 of Tabled Both temperatures show a distinct increase in the region of L183-N1 where the excitation 
temperature is approximately equal to our assumed kinetic temperature. The CLASS hfs fitting assumes an 
equal excitation temperature for all hyperfine components, however, radiat ive transfer rnodeling has shown 



that the temperatures of the separate components can vary by up to 15% dCaselli et al.lll995l : IPagani et al. 



20071 : Daniel et al.ll2007h . This effect is negligible for our purposes as we are primairly interested in the 



kinematics and we are in a low optical depth regime. 
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4.5. Velocity gradients 



Figure [5] shows velocity gradient vectors. These vectors are coherent across the northern half of the 
core - a pattern that can be interpreted as rotation. To investigate this further we compute the the mean 
velocity gradient Q over a 5 by 5 pixel box and the direction of the mean gradient for each of the L183-N 
sources. These parameters and their standard deviations are listed in columns 11 and 12 of Table [2l The 
vectors plotted in Figure |5] have a cut off of ^ > 2 kms^^ pc^^, but all non-zero values were used for the 
calculation of Q. 

Both L183-N1 and L183-N3 have similar magnitudes and directions of velocity gradient vectors. The 
large standard deviation for the magnitude of the L183-N3 gradient is due to the large scatter in the gradient 
vectors, whereas for L183-N1 the large deviation is in part due to an increase in the gradient from the east 
to the west side of the core. The motions across L183-N2 show no ordered direction and the mean position 
angle is not statistically significant. 

Caselli et al. ( 2002 ) and Crapsi et al. ( 2005al') use sing l e dish N2H+ observations to reconstruct the 



Caselh et al 



large-scale velocity field across the L183 core, 
kms^^ pc^^ at a position angle of — 55±3°, while lCrapsi et al 
at a position angle of — 49±1°. 



( 20021) fi tted a velocity gradient of 1.19±0.08 



(l2005al) fitted a gradient of 1 .4±0. 1 km s"^ pc" 



Pagani et al.l (|2007l ) used the IRAM 30m telescope to measure N2H"'" emission at 12-arcsec intervals in 



an EW strip across the filament. They noted a small anti-symmetric velocity shift of a few tens of m s^^ at 
distances further than ±30 arcsec from the center of the core. This gives a large-scale velocity gradient, for 
the cloud as a whole, of 1.2 km s~^ pc^^, which is comparable to the other single dish results quoted above. 

The direction of the large scale velocity field as mapped by the single dish observations is consistent 
with the direction seen in our interferometric measurements. Indeed, the reg ion where we see no gradient 
over L183-N2 is consistent with a region of very small vectors in figure 15 of ICrapsi et al.l (|2005al ). and the 
vectors across the northern part of the core are also in broad agreement. 

While the direction is consistent, the magnitudes of the interferometric and single dish velocity fields 
appear to differ considerably. To check our gradient calculation we can work out the gradient just between 
L183-N2 and L183-N3 using the data in TableU The velocity difference between L183-N2 and L183-N3 is 
0.173 km s~^ and the projected separation between them is 0.018 pc. This gives a velocity gradient between 
them of ' 



' 10 km s ^ pc ^ which is consistent with the 



■^8.9 km s ^ pc ^ measured around L183-N3. 



It appears that the single dish observations are sensitive to a large-scale, low-magnitude velocity gradi- 
ent while the interferometric studies are sensitive to a smaller-scale, larger-magnitude velocity gradient. Yet 
the direction of the gradients on the small scales is consistent with the large-scale gradients. 



The mean direction of the magnetic field in the plane of the sky observed by ICrutcher et all (|2004h 
is 28° ± 2°, whi ch is roughly orthog onal (85° - 102° ) to the velocity gradient vectors in Figure [5] and as 
observed by both Icaselli et all (|2002l ) and lcrapsi et all (|2005al) . 
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4.6. Comparison with single dish observations 



The left hand panel of Figure |7] shows the 5.6 arcsec pixel data smoothed spatially to the resolution 
of a 30 m single antenna telescope. Our BIMA FWHM at 3.22 mm was 13 x 10 arcsec, the FWHM of a 
30 m telescope (e.g. the IRAM dish) at the same wavelength is 27 arcsec. We use the term "30m-Beam" 
to note when we are referring to the smoothed equivalent of a 30 m single-dish beam. The smoothed map 
was multiplied by a factor a 5.6, the ratio of the 30m-Beam to BIMA beam integrals, in order to give a map 
calibrated in Jy/30m-Beam. The conversion factor between Jy/30m-Beam and K becomes 0. 19 1 K/(Jy/30m- 
Beam). The smoothed map shows that the multiple peaks we detect with the BIMA interferometer would 
not have been resolved in integrated emission if the observations had been conducted with a 30-m diameter 
single dish telescope. 

While the L183-N1, 2, & 3 substructure may not have been se en in integrat e d inte nsity maps, it could 
theoretically have been resolved in channel maps. However, the IPagani et all (l2007h observations were 
along an west - east s lice t hrough L183 - N2 an d do not comprise a fully sampled map. Furthermore, the 
Dickens et all (|2000l) and Icaselli et all (|2002[) N2H+ observations of LI 83 were taken with th e FCRAO 
which has an even lower resolution (50 arcsec) beam at this wavelength. Nevertheless, Figure 2 of Dickens et al. 
(|2000 ) shows an offset in the distribution of N2H+ in the 2.4-2.6 km s~^ and 2.6-2.8 kms~^ chan nels that 
is comparable to us seeing higher velocity material to the west of the core. ICaselli et all (|2002|) show a 
comparable vector diagram, albeit over a larger area, but they do not show channel maps. 



The (12", 0") position from Figure I of IPagani et all fcoOTH (15*^54^09.3^ -02°52'48") is within the 
same 5.6 arcsec pixel as our L183-N2 position. We repeated the hfs fitting for the smoothed data at that 
location. The results are shown in the right hand panel of Figure |7l The fitted Visr for the smoothed data is 
2.410 lb 0.002 km s~^ which is comparable to that fitted to the unsmoothed data. 

The linewidth /S.y becomes 0.111 ± 0.004 kms^^ which is narrower than the value fitted to the un- 
smoothed data. L183-N2 is at the systemic velocity of L183, is extended spatially, and shows no velocity 
structure. Therefore spatially smoothing the wings of the emission (e.g., channels 2.29 and 2.52km s^^ in 
Figure m would move flux out of the beam whereas smoothing the extended emission in channels near the 
line center (e.g., channels 2.40 and 2.44km s^^ in Figure|2ll would bring more flux into the beam. This differ- 
ential effect will naturally cause the linewidth to narrow. The opacity becomes r = 3.2 it 0.6, which is higher 
than for the unsmoothed data, but is still substantially lower than reported for single-dish observations. 

The line temperature for our smoothed L183-N2 position was fitted as Tmb = 0.4 it 0.1 which agrees 
with our fit to the unsmoothed data at that position. However, inspection o f Figur e [7] s hows that the fit has 
been unable to recover the full extent of the line temperature. iDaniel et al.l (120070 and IPagani et all (|2007h 
re spectively piiblished spectra towards L183 with T4 and values of 1.5-2.0 K. The (12", 0") position 
of IPagani et all \lQQrk) had a peak ~ 1.8 K. Even allowing for the BIMA main beam efficiency (83%, 
LugtenI 1995 ) our line temperatures are only about a quarter of those found in the single dish observations. 



This agrees with the result of our optical depth fitting to the unsmoothed data. A substantial fraction, 
possibly 75%, of the extended N2H+ emission from L183 is being resolved away by the interferometric 
process. It is this process that has revealed the clumps buried within the extended LI 83 core. 



Right Ascension (300D) Riglit Ascension (SOOD) Right Ascension (3000) Right Ascension (300D) 

Fig. 6. — Images of a) line width, b) optical depth, c) main beam temperature, and d) excitation temperature 
across the mapped area. The overlaid contours are the integrated contours from Figure [T] The range of the 
displayed quantity is shown by the bars at the top of each plot. The three labeled stars show the positions 
listed in Table |2l 




15''54°'10^ 09^ 06^ 5 10 15 

Right Ascension (2000) Velocity (km s"^) 



Fig. 7. — The BIMA interferometric N2H+ L183 data smoothed to the equivalent resolution of a 30m single 
dish telescope. The left panel shows the integrated intensity with contours that start and progress in intervals 
of 13.8 Jy/30m-Beam. The three labeled stars show the positions listed in Table |2] The right hand panel 
shows the spectrum and hfs fit for the smoothed data at the L183-N2 position. Details as per Figured 
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5. Discussion 



We have studied the pre-stellar core LI 83 with high-resolution BIMA interferometric observations, and 
have seen that the core breaks up into three separate clumps, which we have labeled L183-N1, 2 & 3. The 
clump L183-N2 is roughly coincident with the densest central part of L183. Furthermore, the velocity of 
L183-N2 is exactly that of the systemic velocity of the whole of L183. We therefore deduce that L183-N2 
is at the center of LI 83, and that the LI 83 core has formed two further fragments, L183-N1 & 3. 

The virial mass for a uniform density, turbulent, rotating molecular cloud threaded by a uniform mag- 
netic field can be written as M^j^ = Mav + Mb + Mg where the three terms on the right hand side represent 
the turbulent and thermal virial mass, the magnetic critical mass, and the kinematic mass of the system from 
Kepler's third law. This can be written as 



(2) 



where AV is the Unewidth in km s~^, r is the radius of the cloud in parsec, B is the magnetic field strength 
in /iG, and Q i s the velocity gradient du e to rotation in kms~ ^ pc~^. Here we have used the expression 
for Mav from iMacLaren et alj (119881) . ICrutcher et all (12004) used the Chandrasekhar-Fermi method to 
estimate a field strength of 5 = 80 /iG towards L183. The approximate radius of L183-N1 and L183-N2 is 
r ~ 0.01 pc and the radius of L183-N3 is r ~ 0.005 pc. Based on these values and the parameters listed in 
Table|2]we estimate M^ir ~ 0.1 M© for L183-N1 and 2 and M^ir ~ 0.03 Mq for L183-N3. 

For a virialised system the virial mass Mmr will be balanced by the actual mass of the fragment Mdump 
plus a surface term Mp due to a confining medium. The material we have resolved away will act as this 
surface pressure term. By equating an expression for the energy in a confining isothermal pressure to the 
standard formula for gravitational potential energy it can be shown that Mp = 316r^\/"-4^io where is 
the number density of the confini ng medium in un its of 10^ cm~^ and Tio is the gas temperature in units 
of 10 K. Using 714 = 100 from iKirk et all (l2005h and our assumed gas temperature of 7K we estimate 
Mp ~ 0.3 Mq for L183-N1 and 2 and Mp ~ 0.07 Mq for L183-N3. In each case the value of Mp is 2-3 
times greater than the value of M^r. Even without estimating the actual mass of these fragments, we can 
see that they are each pressure confined and will probably each collapse to form a protostar. 

Spitzer archival data of L183 was checked for the presence of a you ng stellar object (YSO). No YSO 



candidates were identifie d towards L183 based on the selection crite ria (IHarvey et all 12007 



Rebull et al. 



20071 : iHarvev etal.ll2008h used by the Spitzer c2d (lEvans etal.ll2003h and Gould Belt (lAllen etal.ll2009h 
legacy surveys. There are no 2MASS sources within the LI 83 fields shown in this paper. We therefore 
confirm that L183 and each of the L183-N clumps are pre-stellar in nature and do not contain embedded 
YSOs. 

The large - scale s tructure of LI 83 is elongated approximately north-south at a position angle of ~ —20°. 
Crutcher et all (120041) gave the position angle of the long-axis of the SCUBA emission as ~ —15°. These 
position angles are consistent with the elongation orientation of L183-N2. We have examined the detailed 
velocity structure of the core and found that L183-N1 & 3 each have velocities slightly offset from that of 
L183-N2 - one positive and one negative. Furthermore, there is a velocity gradient in the general direction 
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from L183-N1 to L183-N3 of ~6 - 9 km s ^ pc ^, over a separation of ~0.016 pc (see Table[2l). We refer 
to this hereafter as the small-scale velocity gradient. 

The mean direction of the small-scale velocity gradient across the entire core is — 60 ± 62°. However, 
the mean direction for "high- velocity vectors" (defined as those vectors with Vur greater than 0.1 kms~^ 
away from the systemic velocity, i.e. \Visr — 2.415| > 0.1 kms~^) is —72 it 2°. The high-velocity vec- 
tors have a more coherent direction than vectors with a velocity close to the systemic velocity. The mean 
direction for the subset of high-velocity vectors in L183-N1 is —68 ± 30° (9 vectors), and for L183-N3 is 
—79 ± 5° (5 vectors). These directions are measured over the same areas used in Section[ 



Therefore the axis of rotation for the small-scale gradient is +22 ib 30° for L183-N1 and +11 ± 5° 
for L183-N3, assuming that the axis is orthogonal to the direction of the velocity gradient. We note that the 
axis of rotation given by the mean of all 14 high- velocity vectors in L183-N1 and N3 is +18 + 19°. This i s 
approximately parallel to the direction of the magnetic field across L183 (+28 ± 2°. ICrutcher et al.ll2004t) . 
but is at an angle of ~ 33° to the core long axis. This is indicating that perhaps the magnetic field may be 
playing an important role in the evolution of LI 83. It also indicates that, for LI 83, the axis of rotation shares 
the same offset from the orientation of the cloud as is found i:o r the magnetic fields in other prestellar cores 
Jcrutcher et aLll2004l : llCirk et al.ll20o3 : Iward-Thompson et al.l [2009). 



We therefore interpret our data of LI 83 as a rotating elongated filamentary core, whose rotation axis 
roughly coincides with its magnetic field direction. L183-N2 is on-axis, and therefore essentially station- 
ary. L183-N1 & 3 lie either side of the rotation axis and have similar magnitude line-of-sight velocities in 
opposite directions. The small-scale velocity gradient therefore tells us about the rotation near the center of 
L183. 



On the larger scale of the LI 83 pre-stellar core as a whole, previous single-dish work by ICrapsi et al. 



(l2005aj) has found a core 0.09 x 0.05 pc in extent - see their figure 15. These authors also saw a velocity 
gradient roughly parallel to the gradient we observed, but on a larger spatial scale. They see a gradient of 
~1.4 kms^^ pc~^, over a scale of ~0.05 pc. We refer to this hereafter as the large-scale velocity gradient. 

The ratio of the spatial scales of the large-scale and small-scale velocity gradients is 0.05/0.016 ~ 3. 
The ratio of the velocity gradients themselves is (6 - 9)/1.4 ~ 4 - 6. Given the uncertainties in all of the 
foregoing estimates and calculations, the agreement of these two ratios to within a factor of two is very 
interesting. This is exactly what would be predicted by the law of conservation of angular momentum if 
L183 were collapsing and spinning up as it collapses. 

We therefore hypothesize that all of our data of LI 83 can be explained in terms of a collapsing, frag- 
menting elongated core, that is rotating about the magnetic field direction and sp inning up as it colla pses. 
We note fina l ly tha t our hypothesis is consistent with the single-dish data of both Icrapsi et all ( 2005a ) and 
Pagani et al.l (|2007|), whose large-scale velocity vectors also appear to increase in magnitude away from the 
center of LI 83. Thus our hypothesis is also consistent with previous data of LI 83. 
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6. Conclusions 

We have used BIMA to make deep N2H+ (1-0) maps of the pre-stellar core LI 83. The interferometric 
process resolved away the majority of the extended N2H+ emission and revealed three clumps which we 
have labeled L183-N1, 2 & 3. L183-N2 is at the estabUshed systemic velocity of L183 and is coincident with 
the position of the peak of the dust emission. We detect a velocity gradient across the core. The direction 
of the gradient is consistent with previously published lower-resolution studies, but its magnitude is larger 
than previously seen. We interpret this pattern as the spin up and fragmentation of a rotating and collapsing 
pre-stellar core. 
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